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SUMMARY 


PROBLEM 

The  ejection  by  solar  flares  of  large  numbers  of  protons  (and  other  charged  particles) 
and  their  direction  into  the  polar  ionosphere  by  the  geomagnetic  field,  effectively  lowers 
the  reflection  height  for  vlf  propagation.  A study  of  these  effects  was  made  using  vlf  data 
obtained  in  a high  altitude  investigation  in  Alaska  during  the  period  between  May  1969  and 
November  1970. 


RESULTS 

Analysis  of  the  data  obtained  during  the  Alaska  investigation,  together  with  pre- 
viously published  information  on  the  relation  of  phase  advances  due  to  solar  protons  and 
proton  flux  data  obtained  from  Explorer  41 , allowed  the  development  of  an  empirical 
model  which  was  used  to  reduce  phase  errors  of  the  Omega  navigational  system  to  within 
10  centicycles.  This  made  it  possible  to  reduce  navigational  errors  in  the  Omega  system 
from  greater  than  10  nautical  miles  to  less  than  one  nautical  mile. 


RECOMMENDATIONS 

Additiona:  propagation  paths  in  the  polar  region  should  be  investigated  in  order  to 
determine  the  existence  of  an  extra  enhanced  D-region  corresponding  to  the  auroral  oval. 
Knowledge  of  such  a propagation  structure  would  be  very  important  to  the  prediction  of 
the  behavior  of  propagation  paths  crossing  the  low-altitude  extremities  of  the  polar  cap  and 
would  enable  the  Omega  system  to  warn  navigators  of  the  existence  of  propagation  phase 
errors  and  would  increase  the  navigational  accuracy  of  the  system. 
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INTRODUCTION 


It  has  heen  known  for  many  years  that  some  solar  Hares  eject  large  numbers  of  pro- 
tons (and  other  charged  particles)  and  that  the  geomagnetic  field  directs  these  particles 
into  the  polar  atmosphere.  These  solar  protons  increase  the  ionization  in  the  polar  D-region, 
causing  increased  electron  densities  at  altitudes  from  40  to  80  kilometers.  These  increased 
densities  effectively  lower  the  “reflection  height’’  for  vlf  propagation. 

Assuming  long-path  conditions,  at  the  low  Omega  frequencies  (10.2  and  13.6  kHz) 
single-mode  propagation  is  in  effect  and  the  lowered  reflection  height  results  in  phase  ad- 
vances. Recently,  Westerlund*  observed  that  the  maximum  phase  advance  of  a transpolar 
vlf  signal  was  proportional  to  the  log  of  the  protor  flux  (>  10  MeV). 

In  this  report,  an  empirical  model  is  presented  which  relates  vlf  phase  advances  to 
incident  proton  flux.  This  model  uses  the  vlf  data  obtained  in  a high-latitude  propagation 
study  by  the  Naval  Electronics  Laboratory  Center  (NELC)  in  Alaska  during  the  period  from 
May  1969  to  September  1970.  The  model  also  uses  information  on  phase  advances  due  to 
solar  protons  for  15  events  previously  reported  by  Martin.-  Proton  flux  data  was  obtained 
from  Explorer  41 . The  empirical  model  deduced  from  these  data  car.  be  used  to  reduce 
Omega  phase  errors  to  within  10  centicvcles  which,  in  turn,  reduces  navigational  errors  from 
greater  than  10  nautical  miles  to  less  than  one  nautical  rnile. 

At  the  present  time,  the  Omega  system  does  not  have  the  capability  of  warr  ing 
liavigatois  uf  ^liasc-  errors  caused  by  solar  disturbances  lor  even  of  warning  ihai  un  error- 
<'i using  event  is  l"  pnjgtcvd  Ftlturv  Id  I he  OTjcgjt  tysHi'm  dmU  rudiiJr 

these  extended  eiror-correcting  modeling  capabilities. 


ANALYSIS 


UNDISTURBED  DIURNAL  AND  SEASONAL  PHASE  VARIATIONS 

The  undisturbed  phase  variations  need  to  be  accurately  modeled  because  the  dis- 
turbance phase  advances  are  deter'  .ined  by  the  difference  between  the  measured  phase 
and  the  usual  undisturbed  phase. 

The  major  phase  variations  are  diurnal.  During  daylight  hours,  the  phase  is  advanced 
relative  to  the  nighttime.  These  diurnal  variations  are  solar-caused  because  the  D-region  is 
formed  basically  by  solar  X-ray  and  Lyman-alpna  ionization.  V'LF  radio  signals  are  reflected 
in  the  D-region. 

Initially,  the  monthly  average  phase  (at  each  hour)  was  expected  to  be  close  to  the 
undisturbed  phase.  During  the  summer  and  winter  months,  this  is  almost  the  case.  But. 
during  the  equinoxes,  the  sunlight  conditions  change  rapidly  (relative  to  a one-month 
period)  and  the  monthly  average  may  not  correspond  to  any  real  day  whatsoever.  This  ef- 
fect is  illustrated  in  figure  I . 

The  Omega  system  has  developed  an  extensive  “skywave”  correction  scheme  to  al- 
low accurate  estimates  to  be  made  of  the  expected  phase.  However,  in  this  study,  a simpler 
method  has  been  developed.  The  diurnal  variations  of  each  path  were  modeled  by  using 


1.  Westerlund,  S.,  F.  H,  Reder,  and  C.  Asom,  “Effects  of  Polar  Cap  Absorption  Events  on  VLF  Trans- 
nissions,"  Planet,  Space  Science,  Vol.  17,  p.  1329,  I9A9. 

2.  Martin,  J.  N.,  Omega  Phase  Variations  During  PCA  Events,  NELC  TR  :.  835 , 17  August  1972. 
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Figure  I.  Monthly  average  phase  versus  real-day  phase. 

only  that  portion  of  the  path  which  was  in  sunlight.  Indeed,  a single  model  is  good  for  all 
seasons  (as  shown  in  figure  2),  indicating  that  the  seasonal  effects  on  the  diurnal  variations 
are  almost  completely  due  to  the  varying  amounts  of  the  path  in  sunlight.  The  amount  of 
sunlit  path  (for  the  two  paths  reported  upon)  may  vary  from  zero  to  50  percent  during  the 
winter,  from  50  percent  to  100  percent  during  the  summer,  and  from  zero  to  100  percent 
in  the  spring  and  fall  months. 

The  diurnal  phase  variations  appear  complex  often  with  apparent  1 2-hour  periods. 
These  variations  are  immediately  explained  by  examining  the  diurnal  variations  in  the 
amount  of  path  in  sunlight,  as  shown  in  figure  3.  This  1 2-hour  period  is  seen  on  paths 
which  lie  near  the  North  Pole  and  is  caused  by  alternating  ends  of  the  path  being  in  sunlight. 

The  unexpected  fact  that  the  phase  variations  are  almost  completely  explained  by 
the  fraction  of  the  path  being  in  sunlight  indicates  that  the  undisturbed  D-iegion  seasonal 
chemistry  changes  (increased  summer-loss  rates)^’  **  are  offset  by  some  other  factors.  All 
that  is  determined  in  this  report  is  that  the  quiet-time  vlf  reflection  height  is  a function  only 
of  the  day-night  conditions. 


3.  Montbriand,  L.  E„  and  J.  S.  Belrose,  “Effective  Electron  Loss  Rates  in  the  Lower  D-Region  During  the 
Decay  of  Solar  X-Ray  Events,”  Radio  Science,  Vol.  7,  No.  I,p  133-172,  1972. 

4.  Larsen,  T.  R.,  Disturbances  in  the  High  Latitude  Lower  Ionosphere,  NDRE  Report  62,  Norwegian  De- 


fense Research  Establishment,  March  1973. 

S.  Argo,  P.  E.,  and  I.  J.  Rothmuller,  Effective  Electron  Loss  Rates  in  the  Polar  D-Region  During  Polar  Cap 
Absorption  Events,  NELCTN  2890,  18  February  1975. 
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Figure  3.  Diurnal  variation  versus  path  in  sunlight. 


PC  A PHASE  ADVANCE  MODELS 

Once  the  polar  cap  absorption  (PCA)  disturbance  phase  advances  (A0  = measured 
phase  minus  the  normal  undisturbed  phase)  were  determined,  a regression  fit  to  log-log 
( 10  flux  (>10  MeV))  was  used  to  generate  a first-order  model  to  the  PCA  variations.  The 
factor  of  10  in  the  log-log  term  was  chosen  to  give  a needed  expansion  in  the  flux  axis. 
Another  factor  might  be  equally  suitable. 

Although  the  basic  shape  of  the  calculated  phase-advance  curve  was  similar  to  mea- 
surements, a distinct  diurnal  pattern  was  observed.  This  is  shown  in  figure  4a.  New  regres- 
sion fits,  using  only  total-day  and  total-night  data,  yielded  an  envelope  to  the  phase  varia- 
tions as  shown  in  figure  4b.  The  model  was  constrained  to  give  zero  phase  advance  for 
some  minimum  flux  (0.3  protons  cm““st_*  sec-' ) under  all  diurnal  conditions.  These 
diurnal  variations  are  fitted  as  changes  in  the  slope  of  the  A <j>  versus  log-log  ( 1 Of)  relation- 
ship in  figure  5.  This  fit  is  reasonable  because  the  quiet  model  was  developed  to  describe 
diurnal  variations  where  the  usual  “quiet”  flux  is  approximately  0.3  protons  cm~“st-'sec~' . 
The  two  paths  showed  distinct  seascnal  variations  with  th-"  winter  day  being  much  more 
sensitive  to  changes  in  proton  flux  tnan  the  summer  day. 

The  phase  advances  versus  log-log  (lOf)  for  the  Norway-to- Alaska  path  are  plotted 
for  three  different  conditions  in  figure  5 summer  day  ( 1 00  percent  of  the  path  in  sunlight), 
spring  and  fall  day  (100  percent  of  the  path  in  sunlight),  and  winter  night  (zero  percent  of 
the  path  in  sunlight).  During  the  summer,  the  Norway-to-Alaska  path  is  always  at  least 
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Figure  5.  Diurnal  variations  as  changes  in  slope  of  A0  versus  log-log (l.Of)  relationships. 
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partially  sunlit  while,  during  the  winter,  the  path  is  ,'ever  totally  sunlit.  The  spring-fall 
periods  are  more  complex.  Extrapolating  the  1 00  percent  and  50  percent  sunlit  spring-fall 
data  to  zero  percent  sunlit  conditions,  results  in  a phase  advance  (zero  percent  sunlit)  versus 
proton-llux  curve  t.iat  approximates  the  winter  zero-percent  sunlit  data.  For  this  reason, 
the  three  seasons  were  combined  into  one  model.  For  the  Norway-to-Alaska  path,  a com- 
plete d:umal  model  was  developed  for  the  fall-winter-spring  seasons.  The  Norway  summer 
path  is  always  sunlit  1 00  percent  of  the  time,  and  so  no  estimate  of  the  summer  nighttime 
phases  could  be  made. 

Except  for  summer,  the  New  York-to-Alaska  path  has  conditions  in  which  there  is 
sunlight  ranging  fron;  zero  to  100  percent  for  most  of  the  year.  This  complete  range  of 
variation  enabled  the  development  of  a complete  diurnal  model  having  the  fall-spring  model 
merged  with  the  winter  model  for  this  path.  The  summer-day  ( 100  percent  sunlit)  and  the 
spring-lall  day  ( 1 00  percent  sunlit)  phase  advances  were  best  Fit  by  a linear  relationship  to 
log-log  ( lOf)  for  the  range  of  proton  fluxes  which  were  observed.  However,  the  extremely 
large  event  of  2 November  I960  produced  a nonlinear  phase  advance  versus  log-log  ( 1 0f> 
relationship  for  the  winter  night.  The  summer  and  spring-fall  paths  may  also  respond  in  a 
nonlinear  manner  during  events  larger  than  those  observed  during  the  corresponding  period. 
To  simplify  the  winter  model,  two  straight  lines  were  fitted  to  the  curves.  Diurnal  variations 
are  predicted  by  weighted  averaging  of  the  measured  (zero  percent)  and  deduced  (100  per- 
cent) sunlit  curves  and  are  shown  in  figure  6. 

Table  I contains  the  models  for  the  New  York-to-Alaska  and  Norway-to-Alaska 
paths  for  summer  and  winter  (spring-fall)  conditions  and  diurnal  variations.  These  models 
provide  an  accurate  fit  to  the  measured  phase  advances  for  the  fifteen  available  events. 


LOG  1 OG  (10F) 


Figure  6.  Diurnal  variations  predicted  by  weighted  averaging 
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TABLE  I . PROPAGATION  MODELS. 


Norway-Ataska 

Summer: 

A<t>  [log  log  (100 + 0.3)  [-20.-  0.33  * % sunlit] 
model  erroi,  a = 3.6  centicycles 

Winter: 

A0  = -[log  log  ( 100  + 0.3]  * 1 .05  * % sunlit  - XN1GHT  * 1 1 - 0.01  * % sunlit) 
XN1GHT  = ( [log  log  (100  + 0.3)  * 33  of  loglog(100  < 2 
| [log log (100 -0.1)  * 170 of  loglog(10f)>  2 
model  error,  o = 8.8  centicycles 

New  York-Alaska 
Summer: 

A 0=  [log  log  (100  + 0.3)  * [-20.  --0.20  * % sunlit) 
model  error,  o = 9.96  centicycles 

Winter: 

A 0 = 'Ice  log  (100  + 0.3)  * [-37.  - 0.43  * % sunlit) 
model  eiror  a = 9.0  centicycles 


DATA  DESCRIPTION 

In  this  section,  the  individual  ever.  >.  ed  in  this  study  will  be  examined.  Figures 
7 through  22  include  all  of  these  events  ana  _ itain  the  measured  and  modeled  phase  ad- 
vances. The  25  September  1969  event  has  a double  peak  in  the  proton  flux  (or.e  on  25 
September  and  another  on  28  Sentember.)  Both  show  up  in  the  Omega  phase  d^ta. 

The  2 November  1969  event  was  the  largest  in  the  data  set.  The  proton  flux  (>  10 
MeV)  exceedrj  !(P  ster-lcm”^sec-^ , while  the  next  largest  event  had  a proton  flux  (>  10 
MeV)  of  10“.  Notice  the  differences  in  the  phase  advances  caused  by  path  conditions  (sun- 
light variations).  The  Norway  path  never  has  more  than  25  percent  of  the  path  in  sunlight 
while  the  New  York  path  was  modulated  by  zero  to  100  percent  of  the  path  being  in 
sunlight. 

The  30  December  1969  event  was  very  small.  The  proton  flux  never  exceeded  2 
ster-*  cnf^sec-' , with  phase  advances  in  the  10-to-15  centicycle  regime.  The  Norway  path 
was  undergoing  a strange,  unexplained  phase  retardation  during  the  event  and  so  (he 
modeled  advances  have  the  correct  shape  but  are  shifted  by  ~1 0 centicycles. 

During  the  28  January  1970  event,  the  Norway  to  Alaska  path  had  a variation  in 
sunlight  between  zero  and  25  percent  while  the  New  York  to  Alaska  path  underwent  changes 
of  from  zero  to  100  percent.  The  Norway  path  (which  has  almost  no  diurnal  variations) 
shows  the  triple-peak  structure  in  the  solar  proton  flux  (peaks  at  0000  hours  UT  on  30  Janu- 
ary and  0000  hours  UT  on  i February).  The  New  York  path  has  diurnal  variations  which 
almost  equal  the  proton  flux  effects. 

The  6 March  1970  event  was  the  first  in  the  set  to  have  a geomagnetic  sudden  com- 
mencement (SC)  during  the  phase  advance.  According  to  some  models  of  PCAs,  there 
should  be  an  expansion  toward  the  equator  at  the  time  of  the  SC.  This  would  show  up  as 
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Figure  7.  PCA  events,  25,  26,  27, 28  and  29  September  1969. 


13 


-98 
-1  00 


6 12 
H0CT69 


6 12 
l*GCT6V 


Id 


U U 
1 60CTe»9 


12 


18 


12 


it 


Figure  8.  PCA  events,  14  and  15  October  1969. 
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Figure  9.  PCA  events,  2, 3, 4 and  5 November  1969. 
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Figure  10.  PCA  events,  24  and  25  November  1969. 
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Figure  12.  PCA  events,  30  and  31  December  1969. 


Figure  n.  PCA  events,  28,  29, 30  and  31  January  1970  and  1 February  1970. 
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Figure  14.  PCA  events, 6, 7, 8, 9 and  lOMarch  1970. 
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Figure  15.  PCA  events,  23, 24, 25, 26  and  27  March  1970. 
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Figure  16.  PCA  events,  28, 29,  30  and  31  March  1970  and  1 April  1970. 
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Figure  17.  PC  A events,  5,6  and  7 May  1970. 


Figure  18.  PCA  events,  30  and  31  May  1970  and  1,2  and  3 June  1970. 
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Figure  19.  PCA  events,  13,  14,  15,  16  and  17  June  1970. 
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Figure  20.  PCA  events,  25, 26, 27, 2^  . *A  29  June  1970. 
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Figure  22.  PCA  events,  21,22,  23, 24  and  25  July  1970. 
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an  increase  in  the  phase  advance  for  the  New  York  path  (it  is  not  completely  withir.  *he 
“polar  cap”).  The  SC  occurred  on  8 March  1970  and  no  anomalous  phase  ac  /nnccs  took 
place  at  that  time.  The  model  does  not  predict  the  phase  advance  on  .he  next  uay  (9  March) 
which  is  seen  on  both  paths.  Since  the  Norway  path  lies  completely  within  the  polar  cap, 
thf  cause  is  probably  something  other  than  expansion  of  the  polar  cap. 

The  data  of  23  March  1970  on  the  New  York  path  could  not  be  reduced  because  of 
tne  large  variations  (not  propagation-related)  and  so  the  second  event.  28  March  1970 
which  is  really  a continuation),  is  more  completely  .iiown.  The  late  March  period  is  on 
the  boundary  of  the  winter  and  summer  models.  Although  the  winter  model  (which  was 
used)  well  explains  the  variations,  the  over  predictions  on  the  New  York  path  on  I April 
1970  may  be  due  to  the  onset  of  summer  conditions. 

The  event  of  30  May  1970  was  the  largest  summer  event  in  the  data  set.  The  scatter 
in  the  New  York  data  increases  later  in  the  summer  until  PCA  events  are  almost  completely 
ob«cured.  The  signal  strength  of  the  New  York  path  is  marginally  adequate  under  the  best 
conditions  (winter).  Hence,  during  a summer  event,  the  signal  may  be  lost  causing  rapidly 
fluctuating  phase  data  (as  in  the  21  July  1970  event). 

Each  of  the  events  of  13  June,  25  June,  and  07  July  1970  have  spike  phase  advances 
at  the  onset  of  the  PCA.  These  do  not  correlate  with  flares  or  SC's  and  are  believed  to  be 
relate  ) io  "ossof  signal  strength. 


DISCUSSION 


The  data  and  models  presented  in  the  previous  section  indicate  that,  by  knowing  the 
general  path  conditions  (diurnal  sunlit  patterns  and  the  season)  as  well  as  the  proton  flux 
(energy  greater  than  10  MeV),  the  phase  variations  during  a PCA  event  can  be  accurately 

c-’1  uated. 

To  extend  the  phase-advance  models  to  transpolar  vlf  propagation  in  general,  the 
amount  of  the  path  which  is  affected  (path  length  within  the  polar  cap)  must  be  determined. 
The  structure  of  the  disturbed  D-region  across  the  polar  cap  cannot  be  determined  using 
data  from  two  paths  and,  for  this  reason  in  this  report,  it  is  considered  to  be  homogeneous. 
Ft  ture  studies  of  PCA  events  and  vlf  pha.'e  advances  should  include  measurements  on  more 
than  just  two  paths  because  any  variations  in  the  altitude  structure  could  have  an  impact 
upon  the  expected  phase  advances. 

Martin"  found  that,  generally,  the  Norway-to-Wales  path  has  detectable  phase  ad- 
vances more  than  one  hour  before  those  of  the  New  York-to-Wales  path.  Martin’s  explana- 
tion is  that  the  Norway  to  Wales  path  passes  closer  to  the  geomagnetic  North  Pole.  This  is 
in  agreement  with  the  Hakura  ^ PCA  model  which  states  that  the  D-region  electron  density 
enhancements  start  at  the  pole  and  spread  rapidly  to  the  equatorial  boundary.  Thus,  the 
Hakura  model  predicts  that  the  Norway-to-Wales  path  has  enhanced  fluxes  (phase  advances) 
before  the  more  southerly  New  York-to-Wales  path.  The  model  used  in  the  s*udy,  which  is 
the  subject  of  this  report,  assumes  a homogeneous  proton  flux  caused  enhancement  over 
the  entire  polar  cap  and  cannot,  therefore,  predict  Martin's  observations.  The  observed 
initial  advances,  however,  were  small  (less  than  10  centicycles)  and,  for  this  reason,  do  not 
constitute  a major  problem. 


6.  Hakura,  K.,  “Entry  of  Solar  Cosmic  Rays  Into  the  Polar  Cap  Atmosphere,”  Journal  of  Geophysical 
Research,  Vol.  72,  p.  1461, 1969 
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The  accurate  modeling  of  the  v!f  phase  advances  makes  it  possible  to  determine 
various  aspects  of  the  polar  ionosphere  such  -a  the  lat  tudmal  extent  of  the  polar  cap  and 
the  altitude  structure  of  the  effective  loss  rate.  Both  of  these  factors  are  extremely  impor- 
tant to  the  evaluation  of  the  disturbance  impact  upon  various  ionospherically  dependent 
systems. 

It  should  be  emphasized  that  single-mode  propagation  is  assumed  in  all  of  the  evalua- 
tions of  data  in  this  study.  Previous  phase  advance  models  made  no  such  assumptions.  The 
NHLC  Waveguide  computer  program  (which  has  been  used  successfully  to  model  vlf  multi- 
mode  propagation  in  a highly  accurate  marker)  was  used  to  verify  that,  for  these  paths  and 
frequencies,  the  propagation  was  single-mode  dominated.  These  paths  were  approximately 
5 000  kilometers  in  length  over  areas  of  poor  conductivity  (land).  Over-water  (good  con- 
ductivity, sea-water)  paths,  even  at  the  low  frequencies  used,  will  have  multimode  effects 
for  up  to  10  000  kilometers. 


LATITUDINAL  EXTENT  OF  THE  POLAR  CAP 

In  the  Hakura*1  model  of  PCA  morphology,  the  initial  enhanced  D-region  is  above 
65°  corrected  geomagneii-  :*titude  (cgl).  The  expansion  phase  following  a geomagnetic  SC 
then  pushes  the  enhancen.  ■;  t to  iess  than  60°  cgl.  The  Norway  to  Alaska  path  is  completely 
above  69°  cgl  and  so  lies  almost  completely  within  Hakura’s  polar  cap.  The  New  York-to- 
Alaska  path,  however,  lies  completely  below  63°  cgl  and  completely  above  52°  cgl  as  shown 
in  figure  23.  Thus,  if  the  Hakura  model  (developed  using  hf  riometers)  was  applicable  to 
the  vlf  response  to  PCA  events,  the  New  York-to-Alaska  path  would  remain  unaffected  until 
the  expansion  phase,  at  which  time  more  than  60  percent  of  the  New  York-to-Alaska  path 
would  be  affected. 

The  analysis  and  model  of  the  fifteen  events  included  in  this  study  indicate  that 
such  is  not  the  case.  The  New  York-to-Alaska  path  is  shown  to  be  affected  immediately 
and  for  the  entire  duration  of  the  event. 

Using  the  modeled  variations,  it  is  \ ossible  to  calculate  the  extent  of  the  vlf  polar 
cap.  The  ratio  of  the  normalized  New  York-to-Alaska  variation  to  the  normalized  Norway- 
to-Alaska  variation,  together  with  figure  23,  yields  the  latitude  directly.  The  normalization 
should  take  into  account  the  wavelength,  the  length  of  the  path,  and  the  propagation  param- 
eters (such  as  ground  conductivity).  Through  the  use  of  the  NELC  Waveguide  program  it 
was  determined  that  the  normalization  was  0.95  percent.  This  means  that,  for  a given 
reflection-height  change  of  the  polar  D-region,  the  ratio  of  the  changes  will  be  approximately 
proportional  to  the  ratio  of  the  path  lengths  inside  the  polar  cap. 

Notice  fl  at  /A0j^  = 1 indicates  that  both  paths  are  entirely  inside  the  polar 
cap.  The  polar  cap  is  at  56°  cgl,  approximately,  when  averaged  over  the  diurnal  variation, 
with  an  apparent  midnight  and  noon  value  of  52°  and  >8°  cgl.  respectively.  This  hints  at 
some  slightly  oval  shape  as  indicated  in  figure  2‘  . 

The  Norway-to-Alaska  model  on  the  New  York-to-Alaska  path,  with  a multiplying 
factor  of  0.80  (corresponding  to  a circle  of  latitude  56°  encompassing  the  changes),  yields, 
however,  almost  exactly  the  same  results  in  fitting  the  noisy  New  York  to  Alaska  data 
(5  = 10.7  centicycles  as  opposed  to  a = 10.7  centicycies).  Any  attempts  to  obtain  more  in 
formation  on  the  shape  of  the  vlf  polar  cap  will  need  significantly  better  (less  ncisy)  data 
than  was  used  in  this  study. 

The  indications  are.  therefore,  that,  using  the  proton  flux  and  diurnal  sunlit  varia- 
tions for  a given  path  (as  well  as  path  lengths  above  56°  cgl),  accurate  predictions  of  the 
phase  advances  experienced  by  that  propagation  path  will  be  possible. 


30 


NO.  OF  PATHS  IN  POLAH  CAP  (NEW  YORK  ALASKA) 


DETERMINATION  OF  ALTITUDE  STRUCTURE  OF  EFFECTIVE  LOSS  RATE 


The  modeled  pha  e advances  of  the  Norway-to-Alaska  signal  during  these  PCA  events 
can  he  used  to  deduce  t’ie  £ ltitude  structure  of  the  disturbed  D-region,  effective  steady-state, 
loss  rates  ( ^ ) for  wintei  d-iy  and  night  conditions  as  well  as  for^summer-day  conditions.^ 

The  effective  loss  rate  is  defined  by  the  relationship  Ne“  = Q/tj/  during  steady-state 
conditions  where  Ne  is  the  electron  density  and  0 is  the  production  rate.  Thus,  if  the  alti- 
tude structure  of  both  Ne(h)  and  Q(h)  is  known,  ^/(h)  can  be  determined. 

The  ionization  rate.  0(h),  has  to  be  calculated  using  the  method  described  by  Argo 
and  can  be  approximated  as  being  constant  in  altitude  and  direct ij  proportional  to  flux 
(greater  than  1C  MeV)  as  shown  in  figure  25  (Q  = 720*  flux). 

The  cLaron  density  profile  must  be  obtained  by  a more  circuitous  route.  Using 
standard  propagation  models,^’  ^ the  reflection  height  can  be  related  to  the  measured  phase 
advance  as  shown  in  figure  26  where  single-mode  theory  is  assumed. 

Again  assuming  Wait’s^  model  of  an  exponential  ionosphere,  the  electron  density  at 
each  reflection  height  is  obtainable  (Ne  = 200  exp  ( 1 .5*(70-h))J . Hence,  the  reflection  elec- 
tron density  and  its  altitude,  shown  in  figure  27,  are  determined  for  each  phase  advance.  Be- 
cause the  phase  advance  is  a function  of  the  proton  flux  (see  figure  5),  reflection  electron 
density  and  altitude  are  uniquely  determined  for  a given  proton  flux  (as  is  the  ionization 
0(h)).  Hence,  by  using  various  flux  values.  Ne~(h)/Q(h)  can  be  calculated,  yielding  the 
(Kh)sin  figure  28,  for  various  conditions  (day/night,  seasonal). 

The  loss-rate  values  determined  by  this  method  agree  with  the  monotonically  de- 
creasing-with-increasing-altitude  quiet-time  loss  rates  given  by  Poppoff. * ® This  differs  from 
the  conclusion  of  Moutbriand  and  Belrose-^  that  the  loss  rate  increases  (with  increasing  alti- 
tude; between  60  and  80  kilometers. 

A definite  seasonal  variation  is  indicated  by  the  data.  The  summer  loss  rates  below 
70  kilometers  are  several  times  larger  than  the  winter  loss  rates.  Larsen^  postulates  that  the 
seasonal  affect  would  be  due  to  the  cold  summer  mesophere  allowing  a larger  degree  of 
hydration.  This  then  would  lead  directly  to  larger  values  of  the  electron-ion  dissociative  re- 
combination rate  coefficient,  which  constitutes  part  of  t^. 


7.  Argo,  P.  E.,  Electron  Production  Rate  of  So’.ar  and  Galactic  Cosmic  Rays  in  the  Lower  Ionosphere  D- 
Region,  NELC  TR  1783, 11  August  1971 

8.  Wait,  J.  R.,  and  K.  P.  Spies,  Characteristics  of  the  Earth’s  Ionospheric  Waveguide  for  VLF  Radio  Waves, 

NBSTN  300.  " ~ “ ~ ~~ 

9.  Davies,  K„  Ionospheric  Radio  Propagation,  MBS  Monograph  30,  1 April  1964. 

10.  Poppoff,  1.  G.,  R.  C,  Whitten,  R.  C.  Guntor,  J.  E.  Evans,  and  E.  G.  Jolci,  “Data-Gathering  Methods 
Based  on  Atmospheric  Measurements,”  De'ense  Nuclear  Agency  Reaction  Rate  Handbook,  Report 
1948H,  March  1972. 


\l>  EFFECTIVE  ELECTRON  LOSS  RATE 


Figure  28.  Effective  electron-loss  rate. 


SUMMARY 


The  effects  upon  vlf  navigational  systems  due  to  ionospheric  changes  caused  by  solar  proton 
events  (i.e.  PC'As)  can  be  accurately  estimated  using  only  the  geophysical  conditions  along 
the  path  (amount  of  path  in  sunlight,  path  length  above  56°  corrected  geomagnetic  latitude, 
and  season)  and  the  instantaneous  proton  flux  (as  measurec  bv  a satebite  such  as  SOLRAD 
XI).  These  estimates  can  reduce  navigational  errors  to  within  one  nautical  mile. 

An  important  outcome  of  this  study  is  the  discovery  that  the  “vlf  polar  cap”  extends 
to  =*56°  which  is  seven  degrees  nearer  to  the  equator  than  the  63°  normally  used.  In  fact, 
this  large  polar  cap  is  at  the  previous  “worst-case”  latitude  which  was  thought  to  occur  in 
the  late  stages  of  Hakura’s  PCA  model. 

This  vlf  polar  cap,  which  is  the  polar  region  in  which  vlf  signals  undergo  phase  ad- 
vances caused  by  solar  proton  events,  differs  from  the  polar  cap  which  is  usually  defined  as 
the  region  toward  the  pole  from  the  auroral  oval.  This  larger  measure  of  the  polar  cap  does 
not  imply  that  the  auroral  oval  is  extended  in  the  direction  of  the  equator  from  the  vlf 
polar  cap.  The  large  size  of  the  polar  cap  in  relation  to  vlf  signals  probably  results  from  the 
fact  that  vlf  propagation  is  affected  at  lower  altitudes  than  are  riometers,  and  the  higher 
energy  particles  doing  the  ionizing  have  lower  cutoff  latitudes.  The  fact  that  the  polar  cap 
extends  seven  degrees  greater  than  previous  measurements^  in  the  direction  of  the  equator 
is  extremely  important  to  global  vlf  navigational  systems  such  as  Omega  because  the  area 
affected  during  a PCA  event  is  much  greater  than  just  the  polar  cap  itself. 
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In  fact,  any  location  having  at  least  one  propagation  path  (to  an  Omega  transmitting 
station)  that  passes  within  the  polar  cap,  will  cause  position  errors  if  it  uses  the  transpolar 
path.  Large  portions  of  the  earth’s  surface  receive  only  one  non-transpolar  signal.  These 
regions  clearly  must  be  affected  by  any  PCA  event.  Figures  29  and  30  show  the  affected 
areas  for  the  present  and  future  Omega  configuration,  while  figure  31  indicates  the  before 
and  after  positioning  errors  during  a relatively  severe  event.  Notice  that  the  affected  region 
(the  Atlantic  Ocean)  is  extremely  important  to  all  Department  of  Defense  users.  It  would 
appear  extremely  important,  therefore,  to  implement  the  phase  corrections  as  described  in 
this  report. 

The  two  propagation  paths  used  in  this  study  allow  an  estimate  to  be  made  of  the 
extent  of  the  polar  cap.  More  paths  would  enable  a study  of  the  polar  D-Region  structure 
ind  would  provide  an  answer  to  the  question  of  whether  there  is  an  extra  enhanced  D- 
RCj.'on  corresponding  to  the  auroral  oval.  This  structure  would  be  very  important  to  paths 
crossing  the  low-latitude  extremities  of  the  polar  cap  and  knowledge  of  it  would  lead  to 
increased  navigational  accuracy. 
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Figure  29.  PCA  event  using  Liberian  station. 
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Figure  30.  PCA  event  using  Trinidad  station. 
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Figure  31 . Before  and  after  positioning  errors  during  a severe  event. 
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